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Information Propagation Speed in Mobile and 
Delay Tolerant Networks 

Philippe Jacquet, Bernard Mans and Georgios Rodolakis 

Abstract 

The goal of this paper is to increase our understanding of the fundamental performance limits of mobile and 
Delay Tolerant Networks (DTNs), where end-to-end multi-hop paths may not exist and communication routes may 
only be available through time and mobility. We use analytical tools to derive generic theoretical upper bounds 
for the information propagation speed in large scale mobile and intermittently connected networks. In other words, 
we upper-bound the optimal performance, in terms of delay, that can be achieved using any routing algorithm. 
We then show how our analysis can be applied to specific mobility and graph models to obtain specific analytical 
estimates. In particular, in two-dimensional networks, when nodes move at a maximum speed v and their density 
v is small (the network is sparse and surely disconnected), we prove that the information propagation speed is 
upper bounded by (1 + 0(v 2 ))v in the random way-point model, while it is upper bounded by 0(^/vvv) for other 
mobility models (random walk, Brownian motion). We also present simulations that confirm the validity of the 
bounds in these scenarios. Finally, we generalize our results to one-dimensional and three-dimensional networks. 

I. Introduction 

Recent research has highlighted the necessity and the significance of mobile ad hoc networks where 
end-to-end multi-hop paths may not exist and communication routes may only be available through time 
and mobility. Depending on the context, these networks are commonly referred as Intermittently Connected 
Networks (ICNs) or Delay Tolerant Networks (DTNs). 

While there is a large body of work on understanding the fundamental properties and performance limits 
of wireless networks under the assumption that connectivity must be maintained {e.g., since the seminal 
work by Gupta and Kumar [6]), there are only few results on the properties of intermittently connected 
or delay tolerant networks {e.g., [2], [10], [11], [13]). Most of the effort has been dedicated to the design 
of efficient routing protocols (see [17] for a survey) and comparative simulations, using specific mobility 
models or concrete traces {e.g., [14], [18]). A complete understanding of what one can expect for optimal 
performance {e.g., through theoretical bounds) is still missing for many realistic models. 

In this context, the objective of the paper is to evaluate the maximum speed at which a piece of 
information can propagate in a mobile wireless network. A piece of information is a packet (of small 
size) which can be transmitted almost instantaneously between two nodes in range. If the network is 
connected {i.e., an end-to-end multi-hop path exists) information moves at a rather high speed, which can 
be considered infinite compared to the mobility of the nodes. 

Part of this work will be presented in "Information Propagation Speed in Mobile and Delay Tolerant Networks", P. Jacquet, B. Mans and 
G. Rodolakis, IEEE Infocom, Rio de Janeiro, Brazil, April, 2009. 
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We consider a network made of n nodes moving in a domain of size A (in two dimensions a square 
area), under the unit disk graph model {i.e., nodes are neighbors when their distance is smaller than one). 
In order to study the properties of DTNs that are relevant to the field of applications, we are interested in 
very sparse networks and we are investigating the case where the node density ^ is small. Indeed, most 
applications for DTNs are required to work for sparse mobile ad hoc networks (e.g., [14], [17], [18]), 
where intermittent connectivity is due to node mobility and to limited radio coverage. In these cases, 
the mobile network is almost always disconnected, making information propagation stall as long as the 
node mobility does not allow the information to jump to another connected component. The information 
is either transmitted or carried by a node (requiring a store-carry-and-forward routing model). Thus, a 
"path" is an alternation of packet transmissions and carriages, that connects a source to a destination, 
and is better referred (from now on) as a journey. Informally, our aim is to find the shortest journey (in 
time) that connects any source to any destination in the network domain, in order to derive the overall 
propagation speed. 

In terms of related work on the information propagation speed in wireless networks, the problem has 
been studied in static networks. Zheng [20] showed that there is a constant upper bound on the information 
diffusion rate in large wireless networks. Recently, Xu and Wang [15] proved that there is a unified upper 
bound on the maximum propagation speed in large wireless networks, using unicast or broadcast. The 
article [8] evaluates analytical upper bounds on the packet propagation speed using opportunistic routing. In 
contrast, our main focus here is to evaluate the information propagation speed in mobile and intermittently 
connected networks. 

Taking into account the node mobility, some recent papers have presented initial results on the theoretical 
properties of intermittently connected networks, e.g., [2], [5], [10], [11], [13], [19]. The papers [5], [19] 
analyze the delay of common routing schemes, such as epidemic routing, under the assumption that the 
inter-meeting time between pairs of nodes follows an exponential distribution. The authors of [13] took 
a graph-theoretical approach in order to upper bound the time it takes for disconnected mobile networks 
to become connected through the mobility of the nodes. This work uses an Erdos-Renyi network model, 
where the node connections are done independently of the actual topology of the network. In this paper, 
we will depart from this model in order to integrate the topological nature of the network, for an instance 
of n mobile nodes, first in a square map of size A connected according to the unit disk graph model, and 
then generalized to a map of dimension D. 

In [2], an interesting model of dynamic random geometric graphs (based on a random walk mobility 
model) leads to the first precise asymptotic results on the connectivity and disconnectivity periods of the 
network. Unfortunately, this methodology cannot be extended to evaluate the fastest possible information 
propagation. In [10], [11], Kong and Yeh studied the information dissemination latency in large wireless 
and mobile networks, in constrained i.i.d. mobility and Brownian motion models. They showed that, when 
the network is not percolated (under a critical node density threshold), the latency scales linearly with 
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the Euclidean distance between the sender and the receiver, while the latency scales sub-linearly in the 
super-critical case where the network is percolated. A question that remains to be answered is to find 
precise estimates on the constant upper bounds of the information propagation speed in intermittently 
connected mobile networks. In [7], the authors present an initial analytical upper bound on the achievable 
information propagation speed in an infinite network model. Here, we present the first analytical results 
in the more realistic (and significantly more difficult) case of a large scale but finite mobile network 
model, in order to prove rigorous upper bounds on the maximum achievable information propagation 
speed. Moreover, we derive our theoretical bounds on a more general mobility model than those used in 
the literature, while we also compare our analytical results with simulations. 
More precisely, our main contributions are the following: 

• we present a new probabilistic model of space-time journeys of packets of information in delay 
tolerant networks; 

• we upper bound the optimal performance that can be achieved using any routing algorithm in finite 
two-dimensional mobile networks and we derive theoretical bounds on the information propagation 
speed, depending on the node density and the network mobility; 

• we generalize our results for bounded multi-dimensional networks; 

• we verify the accuracy of our bounds via simulations. 

The rest of the paper is organized as follows. We first analyze in detail the case of two-dimensional 
networks; in Section |ll} we introduce the network and mobility model, we define the information prop- 



agation speed metric, and we discuss our main results and the methodology. In Section III we present 



the detailed analysis and the proof of our theoretical upper bounds. We derive asymptotic estimates for 



the propagation speed in sparse networks in Section IV We then generalize our results in a more general 



model of multi-dimensional networks in Section [V] We illustrate the behavior of the bounds depending on 



the network and mobility parameters (such as the node density and change of direction rate) in Section VI 



We compare the analytical bounds with simulation measurements in Section VII We conclude and propose 
some directions for further research in Section IVIIIl 

II. Model and Overview of Main Results in Two-Dimensional Networks 
A. Mobile Network Model 

In the two-dimensional case, we consider a network of n nodes in a square area of size A = Lx L. The 
nodes are enumerated from 1 to n. In the next section, we will analyze the case where both n, L — > oo 
such that the node density v = -| tends to a (small) constant. 

Initially, the nodes are distributed uniformly at random. Every node follows an i.i.d. random trajectory, 
reflected on the borders of the square like billiard balls. The nodes change direction at Poisson rate r and 
keep a uniform speed between direction changes. The motion direction angles are uniformly distributed 
between and 2n. When r > 0, we have a random walk model; when r — > oo we are on the Brownian 
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limit; when r — > we are on a random way-point-like model, since nodes travel a distance of order L 
before changing direction. 

The billiard model is equivalent to considering an infinite area made of mirror images of the original 
square: a mobile node moves in the original square while its mirror images move in the mirror squares. 
The fact that a node bounces on a border is strictly equivalent to crossing it without bouncing, while 
its mirror image enters the square. With this perspective, the trajectory of a node is equivalent to a free 
random trajectory in the set of mirror images of the original square, while the nodes remain distributed 
uniformly at random. 

We adopt the unit-disk model: two nodes at a distance smaller than one can exchange information. 
The average number of neighbors per node is therefore smaller (or equal) than ir^. In [16], Xue and 
Kumar have shown that if the average number of neighbors is smaller than 0.074 log n, then the network 
is almost surely disconnected when n is large. In order to study the properties of delay tolerant networks 
in the context of their applications, we need to look at sparse networks. Therefore, we assume that the 
number of nodes n tends to infinity at the same rate as the area of the network domain square (so that 
the node density remains constant), and we investigate the case where the node density is small. 

Since we are interested in computing upper bounds on the best possible information propagation, we 
do not consider here the effects of buffering or congestion. Indeed, we assume that a piece of information, 
i.e., a packet of small size can be transmitted instantaneously between two nodes in range. Even under 
these assumptions, we are able to derive finite bounds on the information propagation speed. We note 
that these assumptions do not affect the validity of our upper bounds, since they correspond to an ideal 
scenario with that respect; this allows us to capture the fundamental performance limit of DTNs based 
solely on the network mobility and topology. Moreover, in the case of very sparse mobile networks, the 
previous assumptions do not impact on the accuracy of our results, since information transmission occurs 
much faster than the speed of the mobile nodes. 

B. Information Propagation Speed and Main Results 

Our main result is the evaluation of a generic upper bound of the information propagation speed 
(presented later in Theorem [TJ in this section), which in turn allows us to obtain specific bounds for 
particular models. 

In order to evaluate the fastest possible information propagation, we establish a probabilistic space-time 
model of journeys of packets of information in delay tolerant networks that contains all possible "shortest" 
journeys: the full epidemic broadcast. We call the information, the beacon. Every time a new node is in 
range of a node which carries a copy of the beacon, the latter node transmits another copy of the beacon 
to the new node. In our model, journeys are expressed as space-time trajectories, since store-carry-forward 
routing also implies that we must take into account the time dimension. 



To prove our main theorem in Section III we decompose the packet journeys into independent segments 



and we evaluate the Laplace transform of the journey probability density. From the Laplace transform, 
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we are able to establish an upper bound on the average number of journeys arriving to a point z before 
a time t, where z is a 2D space vector expressing the spatial distance from the source that emitted the 
beacon. More precisely, we are interested to find when the density of journeys becomes almost surely. 
We notice that a zero probability of reaching a given point in space in a given amount of time implies 
an upper bound on the information propagation speed. In order to evaluate a constant bound, we will 
consider the asymptotic case where the distance from the source and the time both tend to infinity. Hence, 
using our approach, we obtain theoretical bounds on the information propagation speed by computing 
the smallest ratio of the distance over the given time, which yields a journey probability of zero. The 
asymptotic approach must be interpreted in the following sense: we evaluate the information propagation 
speed to a distance which is a large multiple of the maximum radio range. In fact, in Section |VII[ we 



will see that the propagation speed quickly converges to a constant value as soon as the distance between 
the source and the destination is simply larger than the radio range. Similarly, in [10], in the case of 
disconnected mobile networks, the authors show that the information propagation latency scales linearly 
with the distance in the same asymptotic setting. 

Therefore, the concept of propagation speed is probabilistic. To express the previous discussion using 
mathematical notations, let us consider that the beacon starts at time t = on a node at coordinate 
z o = (xo,yo). Let us initially consider (for simplicity) a destination node that stays at coordinate zi = 
(xi,yi). Let q u (z ,z 1 ,t) denote the probability that the destination receives the beacon before time t. A 
scalar s > is an upper bound for the propagation speed, if for all s > s , limg^ ^z , z 1; l Zl ~ z °l j = 
when |zi — z | — > oo, with |.| denoting the Euclidean norm. For example, if we prove that q„(z , z 1; t) < 
exp(— a\zi — z | + bt + c), then quantity - is a propagation speed upper bound. 

Using the previously described methodology, we will prove the following main theorem, which expresses 
our generic upper bound on the information propagation speed in terms of different values of the network 
and mobility parameters. 

Theorem 1: For a network in a square area A = L x L, where the number of nodes n — > oo and L — ► oo 
such that the node density v = -| remains constant, an upper bound on the information propagation speed 
is the smallest ratio: 



q 

min { - with 9 

p,e>o 



P 2 v 2 + \t+ * 1 -r 



„ 2 



p \ " " T V i - a*J'>w 



where v is the maximum node speed, r is the node direction change rate, while J () and 7i() are 

/fc>0V2^ (fc!) 2 



modified Bessel functions (see [1]), defined respectively by: Iq(x) = X)k>o(f ) 2fc (W » ant ^ 



E 



' x\2k+l 



k>Q\2> (fc+l)!fcr 

Remark: As we will see, quantities p and 9 correspond to the parameters of the Laplace transform 
of the journey probabilities. Quantity p is expressed as an inverse of distance and quantity 9 is expressed 
as an inverse of time, therefore the ratio - has the dimension of a speed. 
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Since quantities Iq{x) and -h(x) are both greater than 1, the previous expression has meaning when 
^ < \. Above this density threshold, the upper bound for the information propagation speed is infinite. 
Such a behavior is expected, since it is known that there exists a critical density above which the graph is 
fully connected or at least percolates (i.e., there exists a unique infinite connected component with non- 
zero probability) [12]. The infinite component implies an infinite information propagation speed according 
to our definition. The exact value of the critical density is unknown, although there are known bounds 
and numerical estimates [3]. However, in the context of mobile delay tolerant networks, we are interested 
to analyze the sub-critical case. We note that the critical threshold obtained from our analysis is smaller 
than the critical percolation density. 

Theorem 1 gives a concise upper bound on the information propagation speed, which we will illustrate 



in detail in Section VI In order to give a more intuitive understanding of the fundamental performance 
limits of the information propagation speed, we derive the following corollaries expressing the qualitative 
behavior of the upper bound when the node density tends to 0. This case models very sparse mobile 
wireless networks, which as discussed are of special interest in the context of delay tolerant networks. 

Corollary 1: When nodes move at speed v > in a random walk model (with node direction change 
rate r > 0), and when the square length L — > oo, but such that the node density ^ — > 0, the propagation 
speed upper bound is asymptotically equivalent to O(^P^v). 

It is important to notice that the speed diminishes with the square root of the density v. 

A special case corresponds to r = 0, which is a pure billiard model (nodes change direction only when 
they hit the border). 

Corollary 2: When nodes move at speed v > with r = 0, and when L — > oo, but with node density 
v = -| — > 0, the propagation speed upper bound is (1 + 0(u 2 ))v . 

It turns out that the propagation speed upper bound at the limit is v. This is rather surprising because 
we would expect that the propagation speed would tend to zero when v — > 0. 

We note that the above results do not contradict the results of [13], although they can not be directly 
compared, since a unit disk graph cannot be modeled like an Erdos-Renyi graph. Indeed if nodes A and 
B are connected to a same third node C, then both will be connected with a much higher probability than 
the probability we would had if they were in an Erdos-Renyi graph. On the other hand, our analysis in 
fact confirms the results of [10], which imply that the information propagation speed tends to a constant 
and finite value in intermittently connected networks; our results give the first estimates of this finite 
information propagation speed. 

III. Analysis (Proof of Theorem 1) 
A. Methodology and Journey Analysis 

Our analysis is based on a segmentation of journeys between the source and the destination. Formally, a 
journey is a space-time trajectory of the beacon between the source and the destination. In the following, 
we first decompose journeys into segments (i.e., space-time vectors) which model the node trajectories and 
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the beacon transmissions in Section III-B Our aim is to decompose journeys into independent segments, 
therefore a technical difficulty comes from the dependence in the node emissions and movements (for 
instance, the direction of an emission depends on the direction of the node movement). However, we see 



how we can use an independent segment decomposition in Section III-C , in order to upper bound the 
journey probabilities. We then calculate the Laplace transforms of each individual segment, and, making 
use of the journey decomposition, we deduce the Laplace transform of the probability density of each 



journey in Section [Ill-D for a fixed length sequence of segments. Finally, an asymptotic analysis on the 



journey Laplace transform (for large scale networks), based on Poisson generating functions, allows us to 
compute when the journey probability density tends to zero, and consequently evaluate an upper bound 



on the information propagation speed in Section III-E 



We assume that time zero is when the source transmits, and we will check at what time t the beacon is 
emitted at distance smaller than one to the destination at coordinate z = (x,y). The beacon can take many 
journeys in parallel, due to the broadcast nature of radio transmissions, and the fact that the beacon stays 
in the memory of each emitter (and therefore can be emitted several times in the trajectory of a mobile 
node). In a first approach and in order to simplify, we assume that the destination is fixed; however, we 
will later see that the destination motion does not affect our results. 

We will only consider simple journeys, i.e., journeys which never return twice through the same node. 
This restriction does not affect the analysis, since if a journey arrives to the destination at time t, then 
we can extract a simple journey from this journey which arrives at time t too. 

Let C be a simple journey. Let Z(C) be the terminal point. Let T{C) be the time at which the journey 
terminates. Let P(C) be the probability of the journey C. In the following, we consider a journey as a 
discrete event in a continuous set of all possible journeys in space-time, and we convert the probability 
weight -P(C) to a probability density. 

Assuming that there are n nodes in the network, we call p n (z , zi, i) the density of journeys starting 
from z at time 0, and arriving at zi before time t: 

p n (z ,zx,t) = lim — - S2 P(C) . 

r^O TXT — 

\n-Z{C)\<r,T{C)<t 

B. Journey Segmentation 

Let us consider a journey where the beacon is carried by k + 1 nodes £ , £\ y £ 2 , . . . , £k- The node £q is 
the source. Let £\ be the first node that receives the beacon from the source, £ 2 the node that receives the 
beacon from t lf etc. We call the A;-tuple (£i,£ 2 , • • • ,&k) the journey relay sequence N(C). 

Lemma 1: The probability distribution of the journey C only depends on the cardinality \N(C)\. 
Proof: Since node motions are i.i.d., any node in N(C) can be interchanged with any other node. ■ 
Consequently, we can split the journey C into segments (s , si, s 2 , ■ ■ ■ , Sfc), where the segments Sj are 
random space-time vectors, and where Sj is the space-time vector that starts with the event: "the beacon 
is received by i". In the special case of i — 0, the event is the origin of the journey. 
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To compute the probability distribution of the segments, we notice that Sj corresponds to a space-time 
motion trajectory of mobile node li (the trajectory can be possibly zero if the node immediately retransmits 
the beacon), and a space-time vector of the beacon transmission via radio (where the time component 
of a transmission is zero). Therefore, in order to decompose a journey, we define two kinds of segments 
modeling the described situations: 

• emission segments S e (u, v): the node transmits immediately after receiving the beacon; v is the speed 
of the node that just received the beacon, and u is the emission space vector and is such that |u| < 1; 

• move-and-emit segments S m (u, v, w) = M(v, w) + u: M(v, w) is the space-time vector correspond- 
ing to the motion of the node carrying the beacon, where v is the initial vector speed of the node 
when it receives the beacon and w is the final speed of the node just before transmitting the beacon; 
the vector u is the emission space vector which ends the segment. 

With the following lemma, we prove that the vector u which ends the move-and-emit segments can be 
restricted to unitary segments. 

Lemma 2: In a "fastest" journey decomposition {i.e., with respect to an upper bound on the information 
propagation speed), move-end-emit segments 5(u, v, w) can be restricted to unitary emission vectors: 
|u| = 1. 

Proof: First, assume that ii and £ i+1 are not neighbors when li receives the beacon. The earliest time 
at which £ i+1 will receive the beacon from £i is when both become neighbors, i.e., when their distance 
is just equal to 1; therefore, the emit vector is unitary. Conversely, if ii and £ i+ i are already neighbors 
when £i receives the beacon, then £ i+1 can receive the packet immediately after ^ and the segment would 
be an emission segment instead. ■ 
Since we want to check when a beacon can be emitted at distance less than one from the destination, 
we do not include the last emission in our journey definition; therefore, the last segment s k corresponds 
only to the space-time motion trajectory of node £ k (or simply, a motion equal to zero). 

C. Decomposition into Independent Segments 

In this section, our aim is to decompose journeys into independent emission and move-and-emit seg- 
ments. However, there is a dependence in successive node emissions and movements; for example, a node 
moving faster meets more nodes than a slower mobile node; similarly, the probability of a meeting between 
two nodes is in fact proportional to the relative speed between the nodes, hence two nodes that meet are 
more likely to move in (almost) opposite directions; therefore, the direction of an emission depends on 
the direction of the node movement. To overcome these difficulties, we will in fact work with an upper 
bound on the journey probability densities, and we show that this upper bound can be decomposed into 
independent segments. 

Thus, our objective is to compute an upper bound on P(C), the probability density that a journey C 
exists. For a fixed journey relay sequence N(C) of size k, the probability density is a vector in R 3k . 
Based on the journey decomposition, we have the expression P(C) = p(s \si)p(s 1 \s 2 ) ■ ■ ■ p(s k _i\s k )p(s k ), 
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where p(sj|sj + i) is the conditional probability density of segment Sj, given the next segment Sj+i. We 
have the following expressions for the conditional probabilities, for all possible combinations of emission 
and move-and-emit segments: 

. p(S e (ux, vi)|S , e (u 2 , v 2 )) = P(ui)P(vi)jj; 
this is the probability of emission segment S^Ui, vi), when we know the next segment (here, an emission 
segment): P(ui) is the probability density of U! inside the unit disk (emissions are equiprobable in the 
unit disk, hence P(ui) = gi~n ^~ = 2|ui|), P(vj) is the probability that the node moves at speed vi, 
and -jj- is the density of presence of a node on the second segment (to make the emission possible); there 
is no dependence on the parameters u 2 ,v 2 of the second segment, since the node receiving the packet 
re-emits it immediately to one of its neighbors (there is no new meeting); 

. p(S e (ui, vi)|5 , m (u 2) v 2 , w 2 )) = P(ui)P(vi)^, for the same reason; 

. p(iS , m (ui,v 1) w 1 )| 1 S , e (u2,v 2 )) = P(M(vi,wi))P(ui)P(u 2 )P(v 2 )max{0,ui • (wi - v 2 )}^; 
this is the probability of the move-and-emit segment ,S m (ui, v 1; Wj), when we know the next segment 
(here, an emission segment); quantity max{0,ux ■ (wj — v 2 )}-^ is the average rate at which a node 
carrying the beacon on the second segment enters the neighborhood range of the previous node on the 
radius ui with relative speed v 2 — wi (see Appendix |A|); quantities P(M(vi, wi)), P(ui), P(u 2 ) and 
P(v 2 ) correspond to the probabilities of the respective space and speed vectors (we note the dependence 
on the parameters u 2 , v 2 expressing the probability of the second segment, since the first segment includes 
a node motion, during which the packet is carried before being transmitted to a new neighbor); 

. p(S m (u h vi,Wi)\S m (u 2 , v 2 , w 2 )) = P(M(vi, Wi))P(ui)P(u 2 )P(M(v 2 , w 2 )) max{0, u r (w 1 -v 2 )}^, 
for the same reason. 

From the above, we notice that a journey cannot be directly decomposed into independent segments, 
because of the conditional probabilities. However, recall that in order to derive an upper bound on the 
information propagation speed, it suffices to compute when the probability of a journey becomes zero. 
Therefore, we can instead use an upper bound on the journey probabilities, and check when this upper 
bound becomes zero. Based on the previous expressions, we can upper bound the conditional probabilities. 
Hence, an upper bound of the density of C is P(C), with P(C) = P(si)P(s 2 ) • • ■ P(s k _ 1 )P(sk), and: 

. P(5 e (u,v))=P(u)P(v)i, 

. P(S m (u,v,w)) = P(u)P(M(v,w))2max{i;}± 
where max{v } denotes the maximum node speed. 

Looking at all the above equations, we observe that P(sj) > P(si\s i+ i) for all i and any combination 
of segments. Using the new segments probabilities, we have an upper bound journey model that can be 
decomposed into independent emission and move-and-emit segments. 

D. Journey Laplace Transform 

Let a = ((, 9) be a is an inverse space-time vector: £ is a space vector with components expressed in 
inverse distance units, and 9 is a scalar in inverse time units. We define gfc(cx) as the Laplace transform of 
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the upper bound density of a journey C given that N(C) is fixed of size k. In other words, we have by the 
Laplace transform definition: gk(cr) = E(exp(—a ■ (s + s\ + ■ ■ ■ + Sfc))), under the probability weight P. 
Notice that the exponent is the dot product of two vectors, and that this product is a pure scalar without 
dimension, since (£, 9) is an inverse space-time vector. 

Lemma 3: The Laplace transform of the upper bound journey density given that the relay sequence 
N(C) is fixed and is of length k, satisfies: 

= (2max{ U }^( e --( M ( v - w ) +u ))+ J B(e-^)) fe Q^ 

XjB ( e --M(v,w)^ 

where max{t>} denotes the maximum node speed. 

Proof: This is a direct consequence of the independence of segments in the upper bound journey 
density model, which implies that the journey Laplace transform can be expressed as a product of the 
individual segment Laplace transforms. The first line expresses the Laplace transform of a sequence 
of k emission or move-and-emit segments, while the last term corresponds to the last segment which 
corresponds only to a space-time motion trajectory. ■ 
Let G n {o) be the Laplace transform of the upper bound density of all journeys in a network of size 
n in a square map of area size A. Now, the remaining difficulty comes from the fact that N(C) is not 
known or fixed. To tackle this problem, we define the Poisson generating function: 

yn 

Lemma 4: The following identity holds: 

G(Z,a) = J2h(v)Z k ■ 

k>0 

Proof: This is a formal identity. Quantity gk(cr) depends only on the actual length of the relay 
sequence and not on the nodes that are actually in the relay sequence (from Lemma [j}, thus the Laplace 
transform of the journeys that are made of k segments is ^"^ (^.(cr), since t^ct is the number of distinct 
relay sequences of size k. 

This means that \Z n \ G(Z,a) = J2 k >n (^ky.9k{°) , and 

G(Z,a)e z = E„>o, ft >„P)!*Wf 

U 

Corollary 3: We have (from Lemmas [3] and [4]): 

fi (7 , = E{e -«-M( V , w)) 

<J A Z ') (T J 1 fn r..1 rrv »./l/lv»,U„» , rrV A„\\Z" 

A 



1 - (2max{?)}E(e- ff -( M ( v . w )+ u )) + E(e~(- U )) 
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Therefore, we can now evaluate the Poisson generating function G(Z, a) of the journey Laplace 
transform by combining the segments Laplace transforms. In the following lemma, we evaluate the 
expressions for the Laplace transforms. 

Lemma 5: We have: 

• E(e~^' u ) = 27rWo(|Cl) when u is unitary and uniform on the unit circle, with density u; 

• E(e~^' u ) = Vt^IiQQ) when u is uniform in the unit disk, with density v\ 

• When all speeds are of modulus equal to v, we have _E(e~ CT ' A/ ( v,w )) = , 1 = — ; 

yf (9+t) 2 ~\(\ 2 V 2 -t 

where Ji() and I () are modified Bessel functions. 

Proof: See Appendix |Bj ■ 

E. Information Propagation Speed Analysis 

Our aim is to obtain an estimate of p n (z , zi, i), i.e., the upper bound on the density of journeys that 
start at z at time and end at zi at time t. Let p(Z, Zq, zi,£) be the Poisson generating function of 
p n (z ,zi,t), that is: p(Z, z , zi, t, z) = J2 n Pn( z o, z i^)|r e ~ Z - 

Lemma 6: The generating function p(Z, zo,Zi,t) has positive coefficients. 
Proof: From Lemma |4} 

p(Z,z ,z 1 ,t) = ^2p k (z ,z u t)Z k . 

m 

Hence, we can use the following depoissonization Lemma. 
Lemma 7: When n — > oo: 

p n (z ,zi,t) < z ,zi,t)(l + o(l)). 



Proof: See Appendix [C] ■ 
ij Space-time Asymptotic Analysis: We now evaluate the asymptotic behavior of the journey den- 
sity p v {z-Q, zi, t). With a slight change of notation, we have substituted the node density instead of the 
number of nodes in the network, since in fact we are interested in the limit where n tends to infinity, 
while v = ^ remains constant. From Lemma |7J we see that we can equivalently evaluate the asymptotic 
behavior of the Poisson generating function coefficient p(v, z , z 1; t) (where the number of nodes n tends 
to infinity). 

As we can observe by substituting the expressions of Lemma [5] in Corollary 3 (again with v — ^, 
the asymptotic coefficient (corresponding to the journey density when n — > oo) of the Poisson generating 
function G(u, a), with a = ((,0) a space-time vector, has a denominator K(\£\,6), such that (with 

P = \c\y. 

K(p, 9) = (1 - ~n-h(p)) (V(r + 6) 2 -p 2 v 2 - r) - fovl (p). 
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The key of the analysis is the set /C of pairs (p, 9) such that K(p, 9) = 0, called the Kernel. In fact any 
element of the Kernel {i.e., a singularity of the Laplace transform) can be used to obtain an asymptotic 
estimate of the journey probability density. We denote (p , 9 ) the element of the Kernel that attains the 
minimum value - p . Notice that (p ,9 ) is a function of v = j^. We prove the following lemma. 

Lemma 8: Let u be fixed and 9 1 > 9 . There exists an A 1 such that, when |z| = \zi — z | and t both 
tend to infinity: 

p(v, z , zi, t) < Ax exp(-p |z| + 9it) . 



Proof: See Appendix |Dj ■ 
2) Information Propagation Speed: Let z , and zi be fixed. Let g„(z , z 1; t) be the probability that 
there exists a journey that arrives at distance less than 1 to a destination node at zi before time t. 
Lemma 9: We have the upper bound: 

g„(z ,zi,t)< / p(u,z',z 1 ,t)dz' . 

</|zi-z'|<l 

Proof: By the definition of q u (z , zi, t). ■ 

Therefore, from Lemmas [7J [8] and [9| when L — > oo with z , zi fixed and z = zi — z , we have the 
estimate q u (z ,z 1 ,t) = 0(exp(9it — |z|p )) for all 9\ > 9 . 

Clearly, q u (z , z 1; t) vanishes very quickly when t is smaller than the value such that — p |z| + 9 t = 0, 
i.e., when ^ = -y. This ratio gives the upper bound for the propagation speed. In other words, point 
(p , 6*o) achieves the lowest ratio - in the kernel set /C. By expressing the kernel set K, using the function 
K(p,9) from the previous section, we obtain Theorem 1. 

Remark: We note that this result corresponds to the situation where all nodes speeds are of modulus 
v (as assumed in Lemma [5]). Even if the speeds follow a different distribution, our analysis still applies, 
with the only change occurring in the Laplace transform of the motion vectors (but then the final form 
of Theorem 1 would be different). However, for an upper bound on the propagation speed, it suffices to 
consider v as the maximum node speed. 

To formally complete the proof, we need to address two remaining details: the contribution of the mirror 
images of the nodes (i.e., to account for the nodes bouncing on the borders) and the destination's motion. 
We note here that all node mirror images induce a contribution factor of order exp(6 l it — |z|p — x ), where 
x is the distance of the node from the border of the square network domain (see Appendix [E]); for almost 
all nodes, x is of the order of L, i.e., the edge length of the square, which tends to infinity; therefore 
the contribution of the mirror images is negligible, since it decays exponentially in exp(— L). For the 
destination's motion, it suffices to multiply the journey Laplace transform with the Laplace transform of 
the destination node excursion from its original position, to compute an upper bound on the propagation 
speed. Similarly, the destination's motion also induces a negligible factor (see Appendix [F]). 
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IV. Sparse Two-Dimensional Networks 

A. The Random Walk Model 

Corollary 1: When nodes move at speed v > in a random walk model (with node direction change 
rate r > 0), and when the square length L — > oo but such that the node density -| — > 0, the propagation 
speed upper bound is asymptotically equivalent to 0(y/~^v). 

Proof: Let (p, 9) be an element of the set /C. We have 9 = a/(t + vH(p)) 2 + p 2 v 2 — t, with 
H(p) = yffi y We have: H{p) = ^ + 0(p 2 ), where i/ = % Therefore, 



We obtain the ratio: 



9 = aA 2 + p 2 v 2 - r + J H{p)u + 0{v 2 ). 

\Jt* + p A v A 



6 H(p)u t 

p p T 2 +p 2 v 2 p 



Quantity ^^- + ^r is minimized with value v\J 2uH r ^ - attained at p = V^lEM _ Therefore ^ is minimized 
at value V <J^& + 0(v 3 / 2 ). U 

B. The Billiard Random Way-point Limit 

The billiard limit is equivalent to setting r = 0. 

Corollary 2: When nodes move at speed v > with r = 0, and when L — > oo but with node density 
v = -| — > 0, the propagation speed upper bound is (1 + 0(v 2 ))v . 

Proof: Now, the kernel set fC consists of the points (p,9(p)) where 9(p) = p 2 + Hi(p) 2 with 
Hi (p) = i^K^TJ ) • ^ n m ^ s case ' me u PP er bound speed is proportional to v with a factor of proportionality 



A" p 



equal to ^1 + {^^j where p minimizes Since = ^4vrmin{^} + 0(u 2 ), we get the 

estimate (1 + 0{v 2 ))v , proving Corollary 2. ■ 
These corollaries are useful in order to see more intuitively the behavior of the upper bound of 
Theorem 1, when the node density is small, and consequently to better understand the fundamental 
performance limits of DTNs. Indeed, the case of sparse networks deserves special attention because of 
the potential applications and the necessity to use a delay tolerant architecture. For instance, in the random 
walk model, it is important to notice that the information propagation speed diminishes with the square 
root of the node density v. Furthermore, it is inversely proportional to the square root of the change of 
direction rate of the nodes (changing direction more frequently implies a smaller information propagation 
speed). In fact, the term in the square root in Corollary 1 is proportional to the expected number of 
neighbors that a node meets during a random step. Conversely, in the random way -point model, we notice 
that, surprisingly, the information propagation speed does not tend to with the node density. In this case, 
the upper bound corresponds to the actual maximum speed of the mobile nodes (for instance, halving the 
node speed implies halving the information propagation speed). 
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V. Multi-Dimensional Networks 

In this section, we generalize our bounds on the information propagation speed when the network map 
is in a space of dimension D, from D = 1 to D = 3. This generalizes the case D = 2 treated throughout 
the previous sections. 

The network and mobility model is an extension of the unit disk model described in Section [IIJ Again, 
we consider a network of n nodes in a map of size A = L D , and we analyze the case where both 
n,A — » oo, such that the node density v = tends to a (small) constant. Two nodes at distance smaller 
than 1 can exchange information. Initially, the nodes are distributed uniformly at random. Every node 
follows an i.i.d. random trajectory, reflected on the borders of the network domain like billiard balls. The 
nodes change direction at Poisson rate r and keep a uniform speed between direction changes, while the 
motion direction angles are isotropic. 

The journey decomposition as well as the asymptotic analysis in dimension D = 2 can be directly 
generalized to other dimensions. We note here that the proofs of Lemmas [T] |2l |3j [4j [6j |7J and [9] and 
Corollary [3] hold independently of the network dimension. Therefore, to analyze the propagation speed 
upper bounds in dimension D, we need to adapt the journey Laplace transform expressions (Lemma [5]); 
we generalize the Laplace transforms in the following lemma. 

Lemma 10: For Y D (p, 9), ^d(p) and ^d(p) defined (depending on D) in Table [IJ we have: 

• E(e~^' u ) = uE D (p) when u is unitary and uniform on dimension D, with density v\ 

• ^(e"^'") = v^ D (p) when u is uniform in the unit line, disk, ball (in dimensions D — 1, 2, 3 
respectively); 

• When all speeds are of modulus equal to v , we have _E(e~°"' M(v ' w )) = 1 1 



D 



Yd( P ,8) 



( r+ f))2_ p 2. u 2 



2pv 



T+B+pv 
r-\-B — pv 



*b(p) 



2 cosh(p) 



2 sinh(p)~ 



^(pcosh(p)-sinh(p)) 



TABLE I 

Definition of Y D (p, 6), S D {p) and ^d(p) (depending on D). 



Proof: Equivalently to the proof of Lemma [5J with D— dimensional integration. ■ 
Moreover, we remark that the final result of the asymptotic analysis in Lemma [8] still holds in the case 
of networks in domains of dimension from D = 1 to D = 3. To adapt the proof, it suffices to substitute 



the respective Laplace transform expressions from Lemma 10 (see asymptotic analysis in the appendix in 
Section [D]) and to compute the inverse Laplace transform in space dimension D instead of dimension 2. 
We can thus prove the following more general theorem. 

Theorem 2: In a network of n nodes in a space of dimension D and size A = L D , where n — > oo 
and A — > oo such that the node density v = remains constant, an upper bound of the information 
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propagation speed is the smallest ratio min Pie> o with: 

1 _ 2vuE D (p) 



0. 



Y D (p,9) l-v* D {p) 
where v is the maximum node speed, r the node direction change rate, and the values of Y D (p, 9), Ho(p) 
and ^d(p) are defined (depending on D) in Table |lj 

Remark: From the definition of ^d{), the previous expression has meaning when v < where 
Vd is the "volume" of transmission radius 1: Vd = 2 in 1— D, Vd = tt in 2— D, and Vd = in 3-D. 
Above this density threshold, the upper bound for the information propagation speed is infinite. Such a 
behavior is expected in dimensions 2 and 3, since it is known that there exists a critical density above 
which the network graph percolates, i.e., there exists an infinite connected component. However, a tighter 
analysis in dimension 1 would yield a propagation speed increasing exponentially with the node density, 
in accordance with the size of the largest connected component. 

Proof: Initially, we consider a fixed destination; however, we note that the discussion of the moving 
destination in the appendix (Section [F]) is valid in other dimensions too, therefore the propagation speed 
upper bound remains unchanged if the destination moves as the other nodes. 



Using the new Laplace transforms of Lemma 10 the asymptotic coefficient of the Poisson generating 



function G(u, a) (defined in Corollary |3J, with a = ((, 9) a space-time vector, has a denominator 
K D (\(\,9), such that (with p = |C|): 

js , a\ 1 2vuB. D (p) 



Y D (p,9) l-uV D (pY 
The set K, of pairs (p, 9) such that K D {p, 9) = corresponds to the new Kernel set. Therefore, from the 
new Kernel expression and Lemmas [8] and [9| we obtain the expression for the information propagation 
speed upper bound: the smallest ratio -, with K D (p,9) = 0. 

Again, to complete the proof, we must account for the fact that the nodes bounce on the network 
domain borders, i.e., to add the contributions of the node mirror images as discussed in Section [IlJ in an 
infinite domain of dimension D. According to the analysis of the two-dimensional case (see Appendix [E]), 
the contribution of the mirror images is negligible (in dimension 1 it suffices to consider only the closest 
mirror image, while in dimension 3, we must consider the 4 closest images). ■ 

VI. Slowness of Information Propagation Plots 

To illustrate the behavior of the upper bound for the information propagation speed when the mobile 
density v varies, we define the slowness, i.e., the inverse of the information propagation speed, for which 
our theoretical study now provides lower bounds. Plotting results (obtained by numerical resolution of 
Theorem 2, or, equivalently, of Theorem 1 in the two-dimensional case) of our lower bounds are presented 
in Figures [TJ|2] and [3j in networks of 1, 2 and 3 dimensions respectively, where we consider a unit maximum 
node speed: v = 1ms" 1 . 
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Fig. 1. Theoretical lower bound of slowness versus mobile node density v when r = (left) and r = 0.1 (right), in 1— D networks. 



slowness (x = 0) 




node density V 



^ slowness (x =0.1) 




node density V 



Fig. 2. Theoretical lower bound of slowness versus mobile node density v when r = (left) and r = 0.1 (right), in 2— D networks. 




Fig. 3. Theoretical lower bound of slowness versus mobile node density v when r = (left) and t = 0.1 (right), in 3-D networks. 



Interestingly, in all dimensions, the limit of the information propagation speed when the node density 
tends to zero corresponds to the maximum node speed in the billiard mobility model (r = 0), while the 
propagation slowness is unbounded for small node densities in the random walk model (i.e., when r > 0, 
the information propagation speed diminishes with the node density). 
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We remark that the slowness drops to at v = 1/Vd, with Vd = 2 in 1— D, Vd = 7r in 2— D, and 
Vd = ^ in 3-D: this corresponds to the limit of our model. Recall, that this is a lower bound of the 
slowness (equivalent to the upper bound for the propagation speed). The actual slowness should continue 
to be non-zero beyond v = l/V D . 

Furthermore, in the two-dimensional case (Figure [2]), we notice that the slowness is in 1 — 0{y 2 ) for 
the billiard - random way point limit {i.e., r = 0), confirming Corollary 2; for the random walk, we 
notice that the slowness is unbounded when v — > 0, confirming the 0{^=) theoretical behavior proved in 
Corollary 1. 

VII. Simulations 

In this section, we evaluate the accuracy of our theoretical upper bound in different scenarios by 
comparing it to the average information propagation time obtained by simulating a full epidemic broadcast 
in a two-dimensional network (as described in Section [I]). For all the simulations, we use a unit-disk graph 
model {i.e., a radio range of lm), and the mobile node speed is 1ms -1 . Two commonly used mobility 
models are simulated: the random way-point model (which corresponds to our setting r = as described 



in Subsection IV-B) and the random walk model (which corresponds to our setting r = 0.1). We study 
the two mobility models (Figure [4] for billiard random way-point mobility, and Figure [5] for random walk 
mobility) for different node density v and area values {v = 0.025 on a 80 x 80 square, v = 0.05 on a 
60 x 60 square, v = 0.1 on a 40 x 40 square, respectively). 

In Figures [4] and [5J we depict the simulated propagation time versus the distance (plots), and we compare 
it to the theoretical bound, i.e., a line of fixed slope (in green - bottom). The slope is obtained from the 
analysis in Section III, it represents the slowness illustrated for corresponding density values in Figure [2] 
In the figures, time is measured in seconds, and distance in meters, therefore, the inverse slope of the 
plots provides us with the information propagation speed in ms^ 1 . What is important is the comparison 
of the slopes at infinity. We notice that the measurements very quickly converge to a straight line of fixed 
slope, which implies a fixed information propagation speed. 

Simulations show that the theoretical slope is clearly a lower bound on the slowness, as proved in 
Theorem 1. We also compare the simulation measurements with a second line of fixed slope (red - top). 
This line is provided only for comparison and corresponds to the heuristic situation where we assume 
that node movements and emissions are completely independent (according to the framework of [7] in 
an infinite network). Interestingly enough, the simulations show that the heuristic bound provides an 
accurate slope (the theoretical slope we provide is smaller, since in order to prove a rigorous bound on 
the information propagation speed, we work with an upper bound on the journey probability density for 



the journey decomposition in Section |III-C[ ). 

VIII. Conclusion 

In this paper, we have initiated a characterization of the information propagation speed of Delay 
Tolerant mobile Networks (DTNs) by providing a theoretical upper bound for large scale but finite two- 
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dimensional networks (Theorem [T] and Corollaries 1 and 2) and multi-dimensional networks (Theorem [2]). 
Such theoretical bounds are useful in order to increase our understanding of the fundamental properties 
and performance limits of DTNs, as well as to evaluate and/or optimize the performance of specific routing 
algorithms. The model used in our analytical study is sufficiently general to encapsulate many popular 
mobility models (random way-point, random walk, Brownian motion). We also performed simulations for 
several scenarios to show the validity of our bounds. 

Our methodology and space-time journey analysis provide a general framework for the derivation of 
analytical bounds on the information propagation speed in DTNs. Therefore, future investigations should 
consider extending the analysis to other neighboring models different from unit disk graphs (e.g., quasi- 
disk graphs, probabilistic models), proving tighter bounds (e.g., similar to the heuristic bound we discussed 
in the simulations section), and generalizing to other mobility models or comparing the results with real 
traces. Another interesting direction for further research would be to compare the implications of our 
analysis on the delay of common routing schemes, such as epidemic routing, with the results presented in 
previous work on DTN modeling [5], [19], under the frequently used assumption that the inter-meeting 
time between pairs of nodes follows an exponential distribution. 
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Appendix 

A. Node Meeting Rate 

We consider two nodes moving at speeds Wi and v 2 , respectively. We compute the average rate at 
which the second node enters the neighborhood range of the previous node on the radius Ui (the vector 
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ui is centered on the position of the first node and is of modulus |ui| = 1). The relative speed of the 
nodes is wi — v 2 . The projection of the relative speed on the vector ui equals (ui ■ (wi — v 2 )) Ui. If the 
dot product is positive, the rate at which the second node enters the neighborhood range of the first node 
at Ui, equals U! • (wi — v 2 )^, where quantity \ is the density of presence of the second node. On the 
other hand, if the dot product is negative, the nodes move in such directions that they cannot meet on the 
radius Ui. 

B. Proof of Lemma\5\(Laplace Transform Expressions) 

• The expression E(e~'>' u ), when u is uniform on the unit circle, with density u, is equal to v J Q 27r e'^ cos< ^d(j) = 



In Taylor series, we have utt Sfc( 2 ) 2fc (fc+i)!fcp w hich in turn is equal to z/7Tj|r.Zi(|C|)- 
• We define a carry segment S c (v) as the space-time vector corresponding to a node motion of constant 
speed v, until the node changes direction. Since nodes change direction at Poisson rate r and the speed 
modules is v, we have the expression (with a = ((,9)): E(e~ (T - s ^ ) ) = tY((,6), where Y((,6) = 
r°° e \c\vt-et e -Tt^ _ 1 = (see also [7]). The node motion vector M(v,w) corresponds to 

JO y'( 9+T )2_| C |2 1) 2 ^ 

an arbitrary sequence of carry segments and a final segment, which ends with the reception of the 
information packet by the final destination, instead of a change of direction. Therefore, we have the 
following simple expression inspired from combinatorial analysis: E(e~ a ' M ( v ' wS> ) = jz^y(ce)- This 
is the equivalent of the formal identity = 1 + y + y 2 + y 3 + which represents the Laplace 
transform of an arbitrary sequence of random variables with Laplace transform y (cf. [4]), while 
the term in the numerator corresponds to the final motion vector before emitting to the destination 
(therefore, this term is not multiplied by the direction change rate r). 

C. Proof of Lemma [7| ( Depoissonization ) 
By Cauchy integration (cf. [9]): 

n\ f e z 
p n (z ,zi,t) = — - <b p(Z,z ,z 1 ,t)-^— [ dZ, 

where the integration loop encircles the origin on the complex plane. We take the circle of center and 
radius n: 



2™/ (|C|). 

• The expression E(e~^ '") when u is uniform on the unit disk, with density v, is equal to: 





Therefore, 




= p(n,z ,z 1 ,t)^ : Jo(n). 



Using Stirling and Bessel asymptotics, n! = \j2itnn n e n (l + o(l)) and I(n) 




complete the proof. 
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D. Proof of Lemma [#] (Asymptotic Analysis) 

We prove this lemma for a random mobility model with speed of constant modulus, but generalization 
is straightforward. In this case, with a = (£, 9) a space-time vector: G(v,a) = K ^ e y with K(\C\,9) = 

(1 " ^/iCICD) (v^ + ^-KI 2 ^ 2 - t - 2vrw/ (|C|)) • 

Notice that K ^ ^ is in fact an analytic function of |£| 2 and of 9 with non negative coefficients. One 
must be aware that the quantity \(\ 2 refers to the sum of the square of the coefficients of ( and not to 
the sum of the square of their modulus, and therefore induces an analytical function. Thus, the definition 
domain of G(u,a) contains all tuples a = ((,9) such that (|C|,#) belongs to the definition domain of 



i 



In particular, if there exists a real tuple (pi, #i) in the Kernel /C, then the tuples (p, 9) such that 9?(6>) > 9\ 
and !R(p) < pi belong to the definition domain of K( s p ^ . 

Due to the asymptotic [1] on modified Bessel functions: I a {p) ~ ^/|^cos (| — ^ — | ), we can enlarge 
the definition domain of K ^ p0 - ) to -Re(p) < pi + | log ^(p)], when S(p) is large. 

Since the Laplace transform of space-time density p(u, z, t) is (the denominator comes from 

the fact that we cumulate all journeys arriving at z within t), p(v,z,t) can be expressed via the inverse 
Laplace transform: 

p(u,z,t) = (^[ [ G(u,a)e^ et dC^, 
Am Jdt(e)=9 1 Jsi(0=<i u 

where (Ci,#i) is an y rea l element of the definition domain. 

For our purpose, we will directly take 9\ > 9 and Ci = ~j^} z - Notice that the integration domain in 

£ is an imaginary plane and in 9 an imaginary axis (i.e., integration is 3-dimensional); this explains the 

cubic factor j^s. 

Without loss of generality, we assume that z is co-linear with the first axis: i.e., z = (x, y) and y = 0. We 
denote ( = (( Xl ( y ). The integral in 9 absolutely converges because it is in | 2 , but we cannot conclude the 
same for the absolute convergence of G((, p) which converges only in A- when 3*(C) tends to infinity. To 
this end, we move the integration surface of ( in a suitable way that will lead to an absolute convergence. 
This move is made possible by the multi-dimensional analytical nature of the functions. 

We first define a function p(a) = p + ia for |a| < B for some B > and p(a) = p + ia + | log(|a|) 
for |a| > £>, such that (p(a),6' 1 ) always belong to the definition domain of K } Q , . Second, we define the 
new surface of integration for £ as the the union of Minkowski hyperbolic sections defined by M (a) = 
{£, |(| = p( a ), 3b > : 3?(C) = &Ci}- In other words &(£,) = 0. From the identity C 2 = p(a) 2 + (%(( y )) 2 , 
we get |3fr(C)| > max{K(p(a)), |9f(C)|}. Therefore, 

| / } exp(-z • 0dfl < f M(a) eM-\zM0)\d(\ 

< exp(-5ft(p(o))|z|) 
- \A 

Since exp (-9ft(p(a))|z|) = e- po|z| (a" |z| ), for |a| > 5, the integral converges absolutely in a when we add 
the contribution of all Minkowski sections, as soon as Izl > 2. 
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E. Contribution of Mirror Images 

Let us denote p(u, z, t) as the density distribution of the journeys that connect two points z and Z\ in 
the original square which delimits the network domain, such that zi — z = z, with the journey starting 
at time and arriving at the destination before t. Notice that the journey does not connect to the mirror 
images of zi. To account for the mirror images, when zi = (xi,y%), we need to add the three closest 
images at: z 1 + 2(L — x\)e x , zi + 2(L — yi)e y and zi + 2(L — xi)e :c + 2{L — yi)e y with e x = (1, 0) and 
e y = (0, 1). Adding all possible periodic mirror images, we get the identity: 

p(i/,z ,zi,t) = p{v,z,t) 

+ p(u, z + 2(L - xi)e x , t) + p(u, z + 2(L - yi)^, i) 
+ p(z/, z + 2(L - xi)e z + 2(L - yi)e y , t) 

+ E(j,*)#(o,o) z + 2 ^' Le * + 2/cLe ^ *) 
+ z + 2(L — a;i)e x + 2jLe x + 2kLe y , t) 

+ p{v, z + 2{L - yi)e y + 2jLe x + 2kLe y , t) 

+ p(u, z + 2{L - xi)e x + 2(L - y 1 )e J/ + 2jLe x + 2kLe y , t). 

The dominant terms in the expression of p(u, z , z 1; t) in addition to p(u, z, t) correspond to the three 
closest images. Since we have shown in the previous section that p(v,z,t) decreases exponentially with 
z, from the above identity, the additional factor induced by the mirror images of a given node is of order 
exp(^it — \z\po — x), where x is the distance of the node from the border of the square network domain. 

F. Moving Destination 

We consider that the destination can move as the other nodes, starting at position z at time t — 0. We 
show that the asymptotic propagation speed upper bound does not change when (z, t) tend to infinity. 

For this end, it suffices to multiply the journey Laplace transform with the Laplace transform of the node 
excursion from its original position. The excursion Laplace transform is obtained from the motion Laplace 
transform £'(e _cr ' M(v ' w ' ) ), and it has the expression 1 1 , where Y((, 6) is defined as previously in the 
Laplace transform calculations (see the proof of Lemma [5] in the two-dimensional case, and the definition 



of Yd((,9) in Lemma 10 for the multidimensional case). Therefore, the Poisson generating function 



G*(Z, a) of the new Laplace transform equals G(Z, a) - i 1 (with G(Z,a) the Poisson generating 



function corresponding to a fixed destination, defined in Corollary ph. The function G*(Z,a) has two 



sets of poles, the set JC (described in Section III-E.l) and the new set K! corresponding to the set 



{(p, 9) : 9 = pv — t}, i.e., the roots of the denominator of the excursion Laplace transform. The last set is 
dominated on the right by JC: for all (p, 9') e JC', there is a (p, 9) in JC with 9 > 9' + B and B > 0. Hence, 
the contributions from JC' will be exponentially negligible (of order exp(— Bt)) compared to the main 
contribution from JC, and the propagation speed upper-bound does not change from the value computed 
in Section ITTLEl 
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Fig. 4. Average propagation time (in seconds) versus distance to 
source (in meters), compared with theoretical (green - bottom) and 
heuristic (red - top) slope, for r = and: v = 0.025 simulated in a 
80 x 80 square (top); v = 0.05 in a 60 x 60 square (middle); v = 0.1 
in a 40 x 40 square (bottom). 
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Fig. 5. Average propagation time (in seconds) versus distance to 
source (in meters), compared with theoretical (green - bottom) and 
heuristic (red - top) slope, for r = 0.1 and: v — 0.025 simulated 
in a 80 x 80 square (top); v — 0.05 in a 60 x 60 square (middle); 
v — 0.1 in a 40 x 40 square (bottom). 



